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FULLY AUTOMATED SINGLE-ZONE ELLIPTIC GRID GENERATION 
FOR MARS SCIENCE LABORATORY (MSL) AEROSHELL AND 

CANOPY GEOMETRIES 


Upender K. Kaul 1 2 3 
Ames Research Center 


ABSTRACT 

A procedure for generating smooth uniformly clustered single-zone grids using enhanced elliptic 
grid generation”’ has been demonstrated here for the Mars Science Laboratory (MSL) geometries 
such as aeroshell and canopy, and it has been incorporated in the software. The procedure obviates 
the need for generating multizone grids for such geometries, as reported in the literature. This has 
been possible because the enhanced elliptic grid generator automatically generates clustered grids 
without manual prescription of decay parameters needed with the conventional approach. In fact, 
these decay parameters are calculated as decay functions as part of the solution, and they are not 
constant over a given boundary. Since these decay functions vary over a given boundary, orthogonal 
grids near any arbitrary boundary can be clustered automatically without having to break up the 
boundaries and the corresponding interior domains into various zones for grid generation. 


1. INTRODUCTION 

A smooth and orthogonal grid around arbitrary geometries is invariably generated using grid genera- 
tion techniques based on the solution of partial differential equations. One such powerful technique 
is based on the solution of elliptic partial differential equations. Elliptic grid generation methods 
are generally used to create smooth grids on which accurate numerical solutions (ref. 1) to a given 
physical problem are obtained. This involves the numerical solution of inhomogeneous elliptic par- 
tial differential equations. The inclusion of inhomogeneous terms in these equations allows a grid 
to satisfy clustering and orthogonality properties in the vicinity of specific surfaces in three dimen- 
sions, and in the vicinity of specific lines in two dimensions. Elliptic grids yield more accurate pre- 
dictions, in general, as compared with other grid generation techniques (ref. 1). 


2. FORMULATION 

Two-dimensional (2-D) form of the inhomogeneous elliptic partial differential equations (PDES) 
for grid generation was first proposed by Thompson et al. (ref. 2). These 2-D PDES contained four 
explicit parameters that need to be prescribed by the user. Later, Steger and Sorenson (ref. 3) 
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prescribed a semi-automatic scheme that reduced the requirement for explicit prescription of these 
parameters to the two parameters, called decay parameters. In a subsequent study (ref. 4), the author 
further enhanced this methodology to fully automate the elliptic grid generation process that com- 
pletely eliminated the need for the explicit user prescription of decay parameters. 

Additionally, in the enhanced fully automated methodology (ref. 4), the decay parameters are no 
longer a specified set of constants chosen manually for the four boundaries for 2-D applications, but 
four decay functions, each a function of one independent coordinate variable over a given boundary, 
which are calculated as part of the solution process. This feature makes it possible to cluster a grid 
nonnal to any arbitrarily shaped boundary. 

Below, a three-dimensional (3-D) analog of the earlier 2-D methodology (ref. 2) is given. Then, the 
development of the semi-automatic (ref. 3) and the enhanced fully automated (ref. 4) methodologies 
is briefly discussed. Finally, the extension of the fully automated methodology (ref. 4) for 3-D appli- 
cations is discussed. 

Three-dimensional governing equations for elliptic grid generation are expressed as: 

^xx + ^yy + ^zz = P(£,fi,0 = -afsgn(S - Si) exp-f-filS - Si|}, (1) 

fixx + fiyy + fizz = Q(£,fi,Q = -Ci'sgn(fi - fii) exp{-di|fi - fii|}, (2) 

Cxx + Cyy + Czz = Rfofi.Q = -efsgn(£ - Q exp{-f|^ - Q}, (3) 

where S, fi and £ are generalized curvilinear coordinates, x, y and z are Cartesian coordinates, and 
P(S,fi,Q, Q(S,fi,0, and R(S,fi,0, are inhomogeneous tenns; a;, bi, c,, d;, e; and fi are manually 
selected constants, and the subscript “i” refers to a particular boundary component associated with 
the problem. 

A simplified 2-D form of equations (1-3) is written as (ref. 3) 


^xx + ^yy = -a,'sgn(fi - fii) exp{-di|fi - fii|}, (4) 

fixx + fiyy = -Ci-sgn(fi -fii) exp{-di|fi -fii|}, (5) 


Equations (4) and (5) were used by the authors to semi-automatically generate the 2-D grids with 
appropriate clustering and orthogonality at the walls. However, the decay parameter, d,. needed to be 
set manually, for a given boundary, fi;. 

As mentioned previously, the present fully automatic boundary procedure was proposed (ref. 4) to 
eliminate the need for manual selection of the decay parameters for each boundary. Also, these 
decay functions, no longer constants, but functions of the independent coordinate variables along a 
given boundary, were calculated as part of the solution process. This fully automatic boundary pro- 
cedure has been used successfully for various 2-D complex geometries (ref. 4). 
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In this paper, extension of the formulation for this automatic procedure for 3-D applications is pre- 
sented. Geometries chosen here are the geometries for the Mars Science Laboratory (MSL) aeroshell 
and canopy. The present procedure makes it possible to generate single-zone grids for the aeroshell 
and canopy geometries. 

Using the fully automated approach (ref. 4), for a given boundary, (£>Q, for example, equations 
(1), (2), and (3) are modified here in the context of equations (4) and (5) and are written in the fol- 
lowing fonn. 


Cxx + Cyy + Czz P 3 (Cd"bC) (6) 

where, p 3 (C,P,C) = sgn(C - Q exp{- fi(£,p)IC - Ci|} 

~ (-a 3 ,(q,ri) + a 3 ,i(^,rj) flCdl) (C - CO) sgn(C - Q 
Pxx + Pyy + Pzz = q 3 (CTbC)’ ( 7 ) 


where, q 3 (^,p,C) = -c 3i (^,p) sgn(£ - Q exp{- f^/n^C - 0) 

~ (-c 3 ,i(C/n) + c 3 ,i(C/n) (C - CO) s gn(C - CO 

Cxx+Cyy+Czz = r 3 (C,r|,C), (8) 

where, r 3 (£,p.Q = -e 3i (£,p) sgn(C, - Q exp{- f i (C,'H)lC - UK 

~ (-e 3 ,i(C J 'n) + e 3 ,i(C/n) ffCdl) (C - CO) sgn(£ - Q 

and, where a^, = ak,i(C,P), Ck,i = CkXCdl), and ek,i = Ck.fUp) and where k = 3 corresponds to the 
C, boundary, under consideration. 

Similar expressions hold for the inhomogeneous terms for the c, (k=l) and p (k=2) boundaries. 

The positive decay parameters, b,, di and f,, for the corresponding boundaries, q, p and respec- 
tively, are expressed as parameter functions, bi(p,Q, d;(^,Q, and f,(q,P) in the present approach, and 
the corresponding terms, a 2 /^,Q, c 2 ,i(^,C), e 2 ,i(UC) and ai,i(P,C), ci,i(p,Q, ei,i(p,Q, hold for p and 
c, boundaries respectively. 

Without loss of generality, one can consider the neighborhood of a given ^-boundary segment i, 

C, - C,i > 0. It can be shown that in a selected region on one side of this boundary segment, where 
fi(C,P)(C - CO « 1 , the governing equations and the inhomogeneous terms have the limiting forms as 
given in reference 4. 

Treatment of a boundary segment, - Q < 0, is analogous. Similar PDES hold for regions close to 
^-boundary and p-boundary segments. 


3 



Rewriting, the limiting governing equations near a £ boundary become (ref. 4) 


^xx + ^yy + ^zz - a 3 >i (S,P) fi(^,Tj) (£ - Q = -a 3 >i (S,P) 

fixx + Pyy + Pzz — c 3,i(^,'n) f(^;P) (C — CO = _C 3,i(^>P) 

Cxx + Cyy + Czz - e 3 ,i(tri) f(£,P) (C - CO = -esi^Tl) 

It can be easily seen that the PDES shown previously represent a self-adjoint operator of the form, 

L(0) = div(k grad(0)) - q0 

and, therefore, boundary constraints, valid in the neighborhood of each of the six boundary seg- 
ments, are incorporated by applying Green’s theorem in three dimensions. For example, for the 
£ boundary, the constraint will be given by 


1 s (50/dn) da = j v {-e 3 /£,p) sgn(^ - Q + e 3ji (^,p) f,(q,p) 0} } dx, (9) 

where 0 = (^,p,Q, da is a differential area element, dx is a differential volume element, n refers to a 
direction that is locally normal to a bounding surface S representing a totality of six surfaces includ- 
ing the boundary segments of interest, and V is a volume enclosed by S. This integral-type boundary 
constraint can be used to calculate the decay parameter analog, e 3 /p,Qf(^,p). Similar constraints 
can be used to calculate the decay parameter analogs, ai ; i(r),Qbi(r),Q, and C2,i(T|,C)di(C>C)- When 
expressed in terms of the generalized coordinate, £, the boundary constraint given in equation (9) 
can be written as follows. 

j si da = | s (50/dn) da = j s (dC/dn) da (10) 

The integral 1 si da in equation (10) can be written as an algebraic sum of six integrals, evaluated 
over the indicated boundary segments: 

j si da = 1 ^max I da 1 Tjmax I da + 1 ^max I da 

_ I ^min I da - J iimin I da - J t^min I da (11) 

where the surface configurations qmax, qmin, etc. represent the corresponding boundary segments 
that together make up the surface S. For the first and fourth integral pair, the second and fifth inte- 
gral pair, and the third and sixth integral pair in equation (11), the following respective relations are 
derived. 

j ^1 da = \ 4 (1/Wan) ai 3 [(xrf + yrf + z T1 2 ) (x^ 2 + y^ 2 + z^ 2 )] 1/2 d£ dp. (12) 

I n I da = | n (I/W0C22) OC23 [(x^ 2 + y% + z T 2 ) (x^ 2 + yf + zf)] 1/_ d£, d£, (13) 

1 ^ I da = 1 ^ (1/J) [a 33 (x n 2 + y,f + Zrf) (x^ 2 + y^ 2 + z^ 2 )] 1/2 dp d£, (14) 
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an = J 2 (^x 2 + + ^z 2 ), 

a 2 2 = J 2 (Px 2 + Py 2 + Pz 2 ), 
a 33 = J 2 (^x“ + Cy 2 + ^z). 



0Cl2 J ( ^xfix T ^yHy T ^zfizX 
0C13 = J 2 ( ^xCx + ^yCy + &Q, 

®23 — J ( fix^x T r|y£y + T] z ^ z ), 

where J = J((x,y,z)/(^,r|,Q) is a Jacobian of the transformation (x,y,z) — > &n,Q- 

Equations ( 12)— ( 14) can be used to express the boundary constraints in the computational space 
(generalized variables). The following governing equations in computational space are solved, sub- 
ject to the boundary constraints derived above. 


ail X l; ^ + OC22 Xi rjri + OC33 X i;K + 2{ai2 Xi£r, + an Xj££ + a 23 X; jT |^) = 
-J 2 { p 3 Xi£ + q 3 x i)T1 + r 3 x 1; c }, 
x; = x, y or z. 


A decay parameter analog, such as e 3i i(^,r|)fi(^,r|), may vary with one or more of the generalized 
coordinates, such as (^,T|), rather than being constant; and this variation is determined as part of the 
solution of the grid problem, rather than being prescribed initially by the user. This grid solution can 
be detennined for either a static grid or a dynamically changing grid. Hence, dynamically changing 
grids can be generated automatically without the user intervention. 

The preceding analysis has focused on the neighborhoods of the grid boundary segments. As noted 
in the preceding, in an interior region, far from the grid boundary segments, the defining partial dif- 
ferential equations become homogeneous, and an orthogonal and uniform grid cell distribution is 
obtained which smoothly transitions from the interior to the boundaries. 


3. RESULTS 

Some key elements of enhanced elliptic grid generation methodology are demonstrated through a 
few selected computational grid generation examples. The geometries for planetary entry descent 
and landing systems, such as an aeroshell and a canopy for the Mars Science Laboratory (MSL), are 
considered here. Single-zone grids are generated with the present grid generation procedure which 
are uniformly clustered at the body. With the conventional elliptic grid generation methodology, 
decay parameters are prescribed by the user as constants over a boundary. If the boundary slope is 
discontinuous, then this prescription fails to generate uniformly clustered grids around sharp comers 
and high curvature regions. To avoid this problem in the conventional approach, grids are typically 
decomposed into multiple zones and the grids are generated separately for these zones. 

In the present methodology, decay parameters are automatically calculated as decay functions, as 
part of the solution. Thus, a uniformly clustered grid is generated over an arbitrarily shaped bound- 
ary as a single zone grid. The grid over the MSL aeroshell geometry, shown in figure 1, is generated 
as a single-zone 51x72x61 grid; 72 points in the meridional direction, 51 points in the streamwise 
direction, and 61 points in the direction normal to the body. An axisymmetric cross-section, 51x61, 
of the volume grid around the aeroshell is shown in figure 2. This grid is shown to be clustered at the 
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aeroshell wall. Away from the aeroshell, the grid uniformly stretches to the farfield boundary. No 
clustering requirement is enforced at the farfield boundary. As shown, the grid generated is orthogo- 
nal and uniformly clustered at the aeroshell surface. 

Figure 3 shows the elliptic grid around the aeroshell bounded by a tunnel wall. Therefore, the grid is 
clustered at both the aeroshell and tunnel walls showing uniform clustering at the tunnel wall comer 
points. Conventional methods fail to enforce clustering in the neighborhood of these discontinuities. 
A 3-D slice through the volume grid including the aeroshell surface is shown in figure 4. 

One of the key elements of the enhanced elliptic grid generation procedure is that since the decay 
function varies over the body surface, it can automatically resolve the curvature of the surface in 
accordance with the clustering requirement. It is clearly not possible to prescribe such a decay func- 
tion manually that is required with the conventional elliptic grid generation schemes. 

As an example, a low speed turbulent flow calculation at a Reynolds number of 2xl0 7 was carried 
out with the k-e turbulence model (refs. 5, 6), and the corresponding pressure contours over the 
aeroshell are shown in figure 5. The flow solver (ref. 7) used for this calculation is a low speed flow 
solver. 

Another grid calculation was carried out for the MSL disk-band-gap canopy, and the grid over the 
canopy geometry is shown in figure 6. Again, the grid shown in figure 5 represents a single -zone 
grid. Pressure contours over the canopy in low speed laminar flow are shown in figure 7. 


4. CONCLUDING REMARKS 

A fully automated 3-D elliptic grid generation method has been developed and demonstrated for 
the 3-D MSL geometries such as the aeroshell and canopy. Clustered body-orthogonal grids for the 
MSL aeroshell and the canopy were generated automatically without any user intervention to 
enforce the clustering properties. This makes the 3-D enhanced elliptic grid generator a powerful 
tool for generation of the volume grids around complex geometries, especially the deforming grids 
about the deploying geometries such as MSL canopies and the inflatable aerodynamic decelerators. 
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6. FIGURES 
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Figure 1 . Grid over the MSL aeroshell geometry. 
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Figure 2. Cross-section of the elliptic grid around the MSL aeroshell; outer boundary is farfield. 
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Axisymmetric cross-sectional grid distribution 
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Figure 3. Cross-section of the elliptic grid around the MSL aeroshell and the wind-tunnel wall. 
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Figure 4. A 3-D slice of the grid in the vicinity of the aeroshell. 
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Figure 5. Axisymmetric pressure contours about the aeroshell in low speed 

turbulent flow, Re = 10 7 . 
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Figure 6. A continuous single-zone grid over the MSL disk-band-gap canopy geometry. 
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Figure 7. Pressure contours over the MSL canopy in low speed laminar flow. 
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